#ifndef _ITERATION_H_
#define _ITERATION_H_

#include "Eigen/Dense"
#include "Eigen/Sparse"

Eigen::VectorXd relaxiteration(const Eigen::SparseMatrix<double> &A, const Eigen::VectorXd &b, const Eigen::VectorXd &ini, int times, double omega = 2. / 3);

#endif